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SUMMARY 

An  analysis  is  made  of  the  production  of  sound  by  a  hydrofoil  with  a  Coanda  wall  jet 
circulation  control  (CC-)  device.  Three  principal  sources  are  identified  in  the  vicinity  of 
the  trailing  edge  of  the  hydrofoil.  The  radiation  at  very  low  frequencies  is  dominated 
by  ‘curvature  noise’  generated  by  the  interaction  of  boundary  layer  turbulence  with  the 
rounded  trailing  edge  of  the  CC-hydrofoil;  this  is  similar  in  character  and  magnitude  to  the 
low  frequency  component  of  the  conventional  trailing  edge  noise  produced  by  a  hydrofoil 
of  the  same  chord,  but  with  a  sharp  trailing  edge.  Higher  frequency  sound  is  produced 
principally  at  the  Coanda  jet  slot.  ‘Passive  slot  noise’  is  caused  by  the  ‘scattering’  by  the 
slot  lip  of  near  field  pressure  fluctuations  in  the  turbulent  boundary  layer  of  the  exterior 
mean  flow  past  the  slot.  This  is  of  comparable  intensity  to  high  frequency,  sharp-edged 
trailing  edge  noise.  However,  the  acoustic  spectrum  is  greatly  extended  to  much  higher 
frequencies  if  the  Coanda  jet  is  turbulent;  the  sound  produced  by  the  interaction  of  this 
turbulence  with  the  lip  tends  to  dominate  the  spectrum  at  frequencies  /  (Hz)  greater  than 
about  Uj/h,  where  h  is  the  slot  width  and  Uj  the  Coanda  jet  speed.  Sample  numerical 
results  are  presented  for  a  typical  underwater  application  that  indicate  that  at  this  and 
higher  frequencies  the  slot  noise  can  be  20  dB  or  more  greater  than  conventional  trailing 
edge  noise,  although  the  differences  become  smaller  as  the  thickness  of  the  slot  lip  increases. 
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1.  INTRODUCTION 

A  jet  directed  tangentially  over  a  rounded  solid  surface  tends  to  cling  to  the  surface  in  a 
manner  familiarly  attributed  to  the  Coanda  effect.  The  phenomenon  was  studied  by  Coanda 
[1]  in  the  1930’s,  but  had  been  noted  as  early  as  1800  [2,  3].  The  curved  trajectory  of  a  fluid 
particle  in  the  jet  close  to  the  surface  is  governed  by  a  balance  at  high  Reynolds  numbers 
between  its  centrifugal  inertia  force  and  the  ‘radial’  pressure  gradient,  and  many  such  flows 
are  well  described  analytically  by  methods  of  potential  flow  theory  [4,  5].  However,  viscous 
frictional  forces  at  the  surface  are  necessary  to  produce  the  initial  deflection  of  the  jet 
towards  the  curved  wall. 

A  Coanda  circulation  control  (CC)  device  consisting  of  a  high  speed  two  dimensional 
jet  discharging  from  a  slot  towards  the  trailing  edge  of  an  airfoil  on  the  suction  side  has 
been  extensively  investigated  because  it  can  greatly  increase  the  airfoil  lift  coefficient  [5  - 
12].  When  the  trailing  edge  is  suitably  rounded  the  jet  follows  the  curved  surface  causing 
the  mean  air  flow  over  the  airfoil  to  remain  attached  and  delaying  separation  until  the 
flow  reaches  a  point  on  the  ‘pressure’  side  of  the  airfoil  effectively  corresponding  to  that 
predicted  for  potential  flow  with  the  same  total  circulation.  The  airfoil  can  thereby  attain 
a  maximum  lift  that  is  many  times  larger  than  that  of  a  conventional  airfoil.  In  particular, 
circulation  control  enables  high  lift  to  be  achieved  at  low  flight  speeds  and  at  small  or  zero 
angles  of  attack;  also,  it  can  replace  mechanical  lift  augmentation  systems,  such  as  trailing 
edge  flaps.  A  Coanda  jet  can  deflect  downwards,  around  the  rounded  trailing  edge  of  a  wing, 
the  mean  jet  flow  from  a  turbofan  engine  mounted  above  the  wing  (configured  for  ‘upper 
surface  blowing’)  [5].  This  permits  the  direction  of  thrust  to  be  controlled  by  adjusting 
the  Coanda  jet  speed,  a  procedure  that  allows  almost  instantaneous  ‘thrust  vectoring’  and 
the  attainment  of  the  very  high  lifts  required  for  vertical  flight  without  the  need  to  deploy 
mechanical  high  lift  devices.  Pulsing  the  Coanda  jet  has  also  been  proposed  as  a  means  of 
controlling  airfoil  vibration  [13]. 

But  the  introduction  of  circulation  control  can  involve  large  increases  in  the  noise 
generated  by  turbulence  interacting  with  the  airfoil.  Thus,  upper  surface  blowing  greatly 
increases  turbofan  jet  noise  (by  25  dB  or  more),  caused  by  the  ‘diffraction’  of  jet  shear  layer 
pressures  by  the  wing  trailing  edge;  additional  increases  of  as  much  as  30  dB  are  observed  ’• 
at  high  frequencies  when  CC  is  applied  to  deflect  the  jet  for  thrust  vectoring,  presumably 
produced  by  the  interaction  of  turbulence  in  the  (near  sonic)  Coanda  jet  with  the  jet  slot 
[14].  Carpenter  et  al.  [2,  15]  stressed  the  importance  of  the  jet  slot  as  a  source  of  sound 
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at  supersonic  Coanda  jet  speeds,  for  both  smoothly  contoured  and  ‘stepped’  slots.  For  the 
contoured  slot  the  dependence  of  the  overall  sound  power  on  jet  speed  U  was  approximately 
proportional  to  U 8,  in  accordance  with  Lighthill’s  theory  of  the  generation  of  aerodynamic 
sound  by  ‘free  field’  turbulence  quadrupoles  [16];  for  a  slot  with  a  step  no  clear  velocity 
dependence  was  evident,  but  the  sound  power  at  low  speeds  varied  much  more  slowly  with 
U,  indicating  the  increased  importance  of  surface  interaction  noise  sources  at  these  speeds 

[17). 

This  suggests  that  turbulence-slot  interaction  noise  is  probably  the  most  important 
component  of  the  sound  generated  by  a  Coanda  CC  device  at  low  speeds.  In  underwater 
applications  circulation  control  of  a  hydrofoil  is  a  potentially  attractive  means  of  augmenting 
lift  at  very  low  speeds  and  small  angles  of  attack.  The  analytical  problem  of  predicting 
the  accompanying  acoustic  radiation,  especially  at  higher  frequencies,  is  the  subject  of 
this  report.  At  very  small  Mach  numbers  the  direct  radiation  of  sound  from  turbulence 
quadrupoles  is  usually  of  no  practical  importance.  Our  analysis  identifies  four  distinct 
surface-source  self-noise  mechanisms:  (i)  the  interaction  of  the  turbulent  shear  layers  in  the 
Coanda  jet  with  the  jet  slot;  the  interaction  of  exterior  boundary  layer  turbulence  with 
(ii)  the  jet  slot  and  (iii)  the  rounded  trailing  edge;  and  (iv)  a  weak  surface  interaction  of 
boundary  layer  turbulence  at  the  separation  point  on  the  pressure  side  of  the  hydrofoil. 
Mechanisms  (ii)  and  (iii)  are  found  to  be  important  at  low  frequencies  at  mean  ambient 
flow  speeds  of  a  few  meters  per  second;  higher  frequencies  are  dominated  by  the  Coanda  jet 
source  (i);  separation  noise  (iv)  is  confined  to  very  low  frequencies  and  amplitudes. 

These  conclusions  are  based  on  an  analysis  of  the  production  of  sound  by  a  two- 
dimensional  hydrofoil  with  a  rounded  trailing  edge  fitted  with  a  Coanda  CC  device  close  to 
the  edge.  The  noise  problem  is  formulated  in  §2  of  this  report,  and  a  qualitative  estimate 
is  made  of  the  relative  importance  of  mechanisms  (iii)  and  (iv).  Prediction  formulae  for 
noise  mechanisms  (i),  (ii)  and  (iii)  are  derived  in  Sections  3  and  4.  A  sample  application  of 
these  formulae  is  made  in  §5.  This  includes  a  comparison  with  the  conventional  ‘trailing 
edge  noise’  produced  by  a  nominally  identical  mean  flow  past  the  trailing  edge  of  a  thin 
hydrofoil  of  the  same  chord,  but  with  no  circulation  control  [17].  Of  course,  the  Coanda  jet 
of  a  ‘well  designed’  CC-device  should  be  laminar.  In  that  case  the  jet  would  make  no  direct- 
contribution  to  the  interaction  noise,  although  entrainment  by  the  jet  is  still  likely  modify 
the  amplitude  of  the  sound  generated  by  the  exterior  boundary  layer  turbulence. 
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2.  EQUATIONS  FOR  THE  HYDROACOUSTIC  RADIATION 
2.1  Integral  representations 

Consider  a  two-dimensional  CC-hydrofoil  set  at  zero  angle  attack  in  a  uniform,  high 
Reynolds  number  mean  flow  at  speed  U  in  the  positive  ^-direction  of  the  rectangular 
coordinate  system  (xi,x2,x3)  (see  Figure  1,  where  the  :r3-axis  is  taken  in  the  spanwise 
direction,  out  of  the  plane  of  the  paper).  To  fix  ideas  the  profile  S  of  the  hydrofoil  may 
be  assumed  to  be  nominally  of  elliptic  cross-section,  except  downstream  of  the  spanwise 
Coanda  jet  slot,  where  the  profile  becomes  a  circular  arc  of  radius  The  mean  fluid 
density  and  sound  speed  are  denoted  respectively  by  p0  and  c0,  and  the  mean  flow  Mach 
number  M  =  U/c0  may  be  regarded  as  infinitesimal.  Sound  is  generated  by  the  turbulence 
within  the  boundary  layers  on  S  and  also  in  the  jet  flow  from  the  slot.  At  the  very  low 
Mach  numbers  of  underwater  applications  the  principal  mechanism  of  sound  production  is 
the  diffraction  of  the  near  field  turbulence  pressures  by  variations  in  surface  geometry;  the 
direct  radiation  from  turbulence  quadrupole  sources  (whose  intensity  varies  as  p0U3M5) 
can  be  ignored. 

In  these  circumstances  the  acoustic  pressure  p(x,  t)  at  position  x  and  time  t  can  be 
represented  in  the  form 

f  dG 

p(x,  t)  =  -p0J—  (x,  y,  t  -  r)  •  (ft  A  v)  (y ,  r)d3ydr,  |x|  -»  oo,  (1) 

where  v  and  Cl  =  curl  v  are  respectively  the  fluid  velocity  and  vorticity.  This  formula  is 
applicable  for  low  Mach  number,  high  Reynolds  number  flows,  when  the  contribution  to 
sound  production  by  surface  skin  friction  sources  is  of  no  significance  [17].  G(x,y,t  —  r)  is 
an  acoustic  Green’s  function,  which  satisfies 

(d2/(?0dt2  -  V2)(?  =  S(x  -  y)S(t  -  r),  G(x,y,t-r)  =  0  for  t  <  r,  (2) 

such  that  dG(x,  y \t  —  r)/dyn  =  0,  dG(x,  y;t  —  r)/dxn  =  0  respectively  for  y  and  x  on  S, 
where  yn  and  xn  are  measured  from  S  in  the  direction  of  the  unit  normal  n  (directed  into 
the  fluid). 

Equation  (1)  is  a  useful  representation  of  the  sound  only  when  the  unsteady  velocity 
and  vorticity  are  known,  and  is  applied  below  in  §2.2.  For  predictive  purposes,  however, 
it  is  more  convenient  to  adopt  an  indirect  measure  of  the  turbulent  field  obtained  from 
observations  or  modelling  of  the  unsteady  surface  pressures.  This  can  be  done  by  using  the 
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following  alternative,  but  equivalent,  low  Mach  number  representation  [18] 

p(x,  t)  =  -p0  j£  G(x,  y,  t  -  r)^-( y,  r)d5(y)dr,  (3) 

where  G(x,y,t  -  r)  is  defined  as  in  (1),  and  vn  is  the  normal  component  (directed  into  the 
fluid)  of  the  upwash  velocity  v  of  the  turbulence  for  incompressible  flow,  which  is  formally 
defined  in  terms  of  the  vorticity  by  the  Biot-Savart  formula  [19] 


v(x,t)  =  curl 


r  ft(y,f)cf3y 
Vf  47t|x  -  y|  ’ 


(4) 


the  integration  being  confined  to  the  nonlinear  region  Vs  of  the  flow  outside  the  viscous 
sublayers  on  S  (i.e.  excluding  ‘bound’  vorticity  on  S)  [18]. 

Applications  of  (3)  therefore  depend  on  a  knowledge  of  the  vorticity  Cl.  However,  the 
upwash  velocity  can  also  be  well  approximated  from  measurements  or  by  using  a  suitable 
empirical  model  of  the  surface  pressure  fluctuations  on  S.  Observe  that  v(x,t)  is  the 
‘kinematic’  velocity  induced  by  the  vortex  field  ft  when  the  influence  of  ‘image’  vortices 
within  S  is  ignored. 


2.2  Source  identification 

The  base  level  radiation  from  the  hydrofoil  is  the  sound  produced  when  the  contribution 
from  any  turbulence  in  the  Coanda  circulation  control  jet  is  ignored.  If  turbulence  levels 
in  the  incident  mean  flow  are  negligible  the  primary  sources  are  then  associated  with 
the  growth  of  the  turbulent  boundary  layers  on  S.  The  sound  from  this  source  can  be 
estimated  independently  of  the  contribution  from  the  Coanda  jet,  and  consists  of  sound 
generated  when  the  turbulence  is  swept  past  the  surface  irregularity  formed  by  the  slot, 
by  its  interaction  with  the  ‘trailing  edge’  of  the  hydrofoil  (which  in  this  case  means  the 
region  of  maximum  surface  ‘sharpness’  or  curvature),  and  also  by  the  surface  reaction  force 
produced  by  the  rapid  acceleration  of  turbulence  away  from  the  hydrofoil  in  the  vicinity  of 
the  separation  zone.  Of  these  mechanisms  the  slot  interaction  source  will  occupy  the  higher 
frequency  portion  of  the  acoustic  spectrum,  and  is  considered  in  detail  in  §4.  The  curvature 
and  separation  sources  are  significant  only  at  low  frequencies  ~  U/l ft. 

The  relative  importance  of  these  two  low  frequency  sources  can  be  estimated  by 
considering  their  contributions  to  the  sound  generated  by  a  spanwise  line  vortex,  imagined 
to  be  inserted  into  the  boundary  layer  on  the  suction  side  of  the  hydrofoil  far  upstream  of 
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the  slot  and  then  to  be  convected  (at  the  mean  flow  velocity)  in  the  mean  flow  around  the 
trailing  edge,  as  indicated  in  Figure  2.  For  the  purpose  of  illustration  it  will  be  temporarily 
assumed  (in  this  Section  only)  that  the  frequency  is  sufficiently  small  that  the  hydrofoil 
chord  may  be  regarded  as  acoustically  compact.  This  may  or  may  not  be  true  in  practice, 
but  the  following  discussion  will  still  supply  an  excellent  indication  of  the  relative  strengths 
of  these  sources. 

Let  Xr(t)  =  (XrijXre)  denote  the  position  of  the  vortex  in  the  aqa^-plane,  and  write 

n  =  -n36(x  -  Xr (t)),  X  =  (xi,x2),  |x3|  <  L/2,  (5) 

where  i3  is  a  unit  vector  in  the  spanwise  (x3-)  direction,  and  T  >  0  is  the  vortex  strength 
(corresponding  to  circulation  in  the  clockwise  sense  indicated  in  the  Figure).  At  the  ‘ends’ 
x3  =  ±L/2,  the  vortex  may  be  assumed  to  bend  down  to  the  surface  S  or  into  an  upstream 
direction,  in  such  a  manner  that  end  effects  may  be  ignored. 

When  the  hydrofoil  is  compact  the  dominant  component  of  the  sound  generated  by 
the  vortex  is  equivalent  to  that  produced  by  a  dipole  source  whose  strength  is  equal  to 
the  unsteady  induced  lift,  and  can  be  calculated  from  equation  (1)  by  substituting  the 
representation  (5)  for  Cl  and  by  using  the  compact  approximation  to  the  Green’s  function, 
as  described  in  detail  in  §3.2  of  [17].  This  calculation  yields 

«<“•*)  1x1  ^°°-  (6) 

where  the  term  in  the  square  brackets  is  evaluated  at  the  retarded  position  Xp(t  —  |x|/c0)  of 
the  vortex  ,  0  =  cos_1(x2/|x|)  is  the  angle  between  the  radiation  direction  and  that  of  the 
mean  lift,  and  the  function  T(x)  is  determined  by  the  cross-sectional  shape  of  the  hydrofoil, 
and  represents  physically  the  stream  function  of  a  hypothetical  incompressible,  irrotational 
flow  past  the  hydrofoil  that  has  unit  speed  in  the  x2-direction  at  large  distances  from  the 
hydrofoil.  Evidently,  sound  production  by  the  vortex  will  be  maximal  at  those  locations 
on  its  trajectory  around  the  hydrofoil  where  its  acceleration  across  the  streamlines  of  the 
hypothetical  potential  flow  is  greatest. 

Figure  2b  depicts  a  typical  prediction  of  equation  (6).  In  performing  this  calculation,  the  • 
stream  function  'F  and  the  trajectory  of  the  vortex  have  been  determined  by  taking  the 
hydrofoil  cross-section  to  be  an  ellipse  of  semi-major  and  minor  axes  respectively  equal  to 
a  and  b  with  b/a  =  0.2.  The  general  case  in  which  explicit  account  is  taken  of  the  circular 
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profile  of  the  CC  hydrofoil  at  the  trailing  edge  requires  a  preliminary  numerical  calculation 
to  determine  T(x),  but  the  rounding  has  been  ignored  in  evaluating  (6),  because  predictions 
of  sound  generation  by  the  vortex  are  only  minimally  influenced  by  these  geometrical 
differences  when  the  trajectory  and  'I'  are  both  approximated  in  this  way. 

The  trajectory  of  the  vortex  T  is  shown  in  Figure  2a  as  a  function  otUt/a  when  the 
hydrofoil  is  at  zero  angle  of  attack,  and  the  lift  produced  by  the  Coanda  jet  (not  shown 
in  the  Figure)  is  characterized  by  a  lift  coefficient  Cl  -  3.5.  This  corresponds  to  a  mean 
circulation  about  the  hydrofoil  equal  to  CloXJ ,  and  the  mean  flow  convecting  the  vortex 
around  the  trailing  edge  has  been  calculated  using  the  potential  theory  of  steady  flow  past 
the  ellipse  with  this  circulation  [19].  We  have  taken  U  —  2.5  m/s,  and  the  chord  2 a  to  be  3 
m.  The  vortex  is  released  on  the  suction  side  of  the  hydrofoil  at  the  midchord  position  at 
a  distance  0.0024a  from  the  surface;  this  distance  is  approximately  equal  to  the  boundary 
layer  displacement  thickness  St  when  a  turbulent  boundary  layer  on  S  is  assumed  to  start 
near  the  leading  edge  of  the  hydrofoil.  In  the  Figure,  the  trajectory  has  been  shifted  to  the 
right  by  2.5%  of  the  semi-major  axis  a  in  order  to  make  the  vortex  path  adjacent  to  the 
hydrofoil  visible 

The  nondimensional  acoustic  pressure  signature 

_ P(x,  t) _ 

(, p0t/2)  (rL  cos  0/coa|x|) 

determined  by  equation  (6)  is  plotted  in  Figure  2b  as  a  function  of  the  nondimensional 
retarded  time  U[t]/a,  [<]=£—  |x|/c0  (corresponding  to  retarded  positions  shown  in  Figure 
2a) .  A  comparison  of  the  figures  shows  that  no  sound  is  produced  until  the  vortex  arrives  at 
the  point  of  maximum  curvature  O.  This  ‘curvature  noise’  pulse  has  a  duration  ~  0.1  a/U 
sec.,  corresponding  to  a  minimum  principal  frequency  ~  10 Ufa  ~  16  Hz.  The  vortex 
trajectory  separates  from  the  hydrofoil  close  to  the  stagnation  point  A;  Figure  2b  shows 
that  the  ‘separation  noise’  produced  by  the  vortex  is  of  very  low  frequency  and  amplitude. 
Therefore,  this  contribution  to  the  overall  radiation  is  henceforth  ignored. 
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3.  CURVATURE  NOISE 
3.1  Calculation  of  the  acoustic  pressure 


In  this  section  the  frequency  spectrum  of  the  curvature  noise  generated  by  turbulence 
in  the  boundary  layer  flow  over  the  rounded  edge  O  of  the  hydrofoil  is  determined  from 
equation  (3),  by  expressing  the  upwash  velocity  v  in  terms  of  an  empirical  model  of  the 
local  turbulent  boundary  layer  wall  pressure.  At  high  frequencies  it  is  not  permissible  to 
assume  (as  in  §2.2)  that  the  chord  2 a  is  acoustically  compact.  Since  these  higher  frequency 
components  are  probably  of  more  interest  in  applications,  detailed  calculations  will  be 
reported  only  for  the  case  of  non-compact  chord.  Low  frequency  estimates  can  subsequently 
be  corrected,  where  necessary,  for  the  influence  of  finite  chord,  by  applying  the  interpolation 
formula  given  in  §4.2  of  [20]. 


Take  the  origin  of  coordinates  at  a  point  O  on  the  trailing  edge  of  the  hydrofoil  (see 
Figure  3),  and  consider  the  Fourier  component  p(x,  co)e~tut  of  the  far  field  acoustic  pressure 
p(x,  t).  When  the  chord  is  non-compact  it  is  shown  in  [18]  that  the  frequency  domain  form 
of  equation  (3)  is 


p(x,  w) 


p0ojy/uT0  sin ^  ^>sin(0/2) 


fs$*(y)vn(y,u))eiKolx  hV3ldS(y),  |x|  -4  oo,  (7) 


7r\/2n\x\ 

where  k0  —  u/c0  is  the  acoustic  wavenumber  and  un(y, u>)  =  ^  J^°00vn(y,t)eiutdt.  The 
function  <E>*(y)  is  a  solution  of  Laplace’s  equation  whose  gradient  is  the  velocity  of  a 
potential  flow  around  the  trailing  edge  (in  the  anticlockwise  sense  in  Figure  1),  and  is 
normalized  such  that 


$*(y)  -¥  Re  {-iz*}  for  \z\  »  5ft,  z  =  yi+iy2.  (8) 

In  terms  of  polar  coordinates  ( t/i , y2 )  =  r,(cos0',sin0'),  the  limiting  form  Re  {— iz*}  = 
yfp sin(#'/2);  this  is  the  potential  for  anticlockwise  irrotational  flow  around  the  edge  of  the 
half-plane  yi  <0,  y2  =  0.  The  angles  6,  ip  define  the  far  field  radiation  direction  x/|x|, 
as  indicated  in  Figure  3,  where  (xi,x2)  =  r(cos$,  sin0)  and  ip  is  the  angle  between  the 
spanwise  direction  (the  a^-axis)  and  the  radiation  direction. 

To  evaluate  the  integral  in  (7)  we  put  dS  =  dsdyz,  where  s  is  the  curvilinear  distance 
from  O  measured  on  S  in  the  aqa^-plane,  as  indicated  in  Figure  3.  If  we  also  introduce  the 
Fourier  integral  representation  on  S 

vn(y,w)  =  Jf  vn(k,k3,u)el^s+k3y3')dkdk3,  (9) 
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then,  as  |x|  — >  oo, 

§*{y)vn(y,w)eiKo\x~hy^dS(y)  =  27re,'Co|x|  JJ  vn(k,  k0  cos  ip,  u)$*(s)etksdsdk,  (10) 

where  $*(y)  =  $*(s)  on  S.  The  principal  contribution  to  the  s-integral  is  from  the 
immediate  vicinity  of  s  =  0,  and  it  is  only  the  behavior  of  $*(s)  near  this  point  that 
is  important.  Therefore,  without  loss  of  accuracy  $*(s)  can  be  approximated  by  the 
corresponding  potential  $*(s),  say,  of  flow  past  a  parabolic  edge  whose  minimum  radius  of 
curvature  (at  s  —  0)  is  5?,  i.e.  by 


K(y)  =  \/(yi  +5R/2)2  +  y| 


sin 


(i  tan'1 

y2  ]\ 

V2 

.v'fei  +  3V2)2  +  »iJ/ 

(ii) 


In  a  typical  application  3?  ~  0.08a,  a  =  1.5  m,  and  the  hydrodynamic  wavenumber 
k  ~  u)/U0  where  U0  ~  10  m/sec.  Thus,  kl ft  ~  0.08wa/17o  >  1  when  the  frequency  /  =  cj/2tt 
satisfies 

/  >  ~  13  Hz,  for  a  =  1.5  m. 

0.5a 

It  may  therefore  be  assumed  in  evaluating  the  integration  with  respect  to  s  in  (10)  that 
Mi  »  1  for  all  acoustic  frequencies  of  interest.  In  this  asymptotic  limit  we  find 

iy°(*)°ik’ds*T\[w\e~m>2' 


and  equation  (7)  becomes 

— p0u)\f%Ko  sin^  ip  sin(0/2)e*'c"lxl  r°°  vn(k,  k0  cos  rp,  a;)e-lfclR/2 


p(x,w) 


y/2\z 


L 


*#i 


dk,  |x|  — >  oo.  (12) 


Now  when  /c3 i  »  1,  the  hydrodynamic  length  scale  is  much  smaller  than  the  surface  radius 
of  curvature,  and  we  can  estimate  the  value  of  vn(k,  K0cosxp,co)  in  terms  of  the  boundary 
layer  turbulence  surface  pressure  fluctuations  by  assuming  the  surface  S  to  be  locally  flat 
(on  the  scale  of  1/A:).  The  pressure  on  S  may  be  imagined  to  have  two  components:  the 
first  is  the  ‘incident’  pressure  pi  produced  by  the  turbulence  when  its  interaction  with  S  is 
ignored;  this  satisfies  p0dv/dt  =  —Vpi  on  S,  where  v  is  the  upwash  velocity.  The  second 
is  the  pressure  pn  ‘reflected’  from  S.  On  a  locally  flat  surface  S:  pi  =  Pr,  where  the  two 
pressures  combine  to  form  the  ‘blocked  surface  pressure’  ps(k,  n0  cos  ip,  u),  say.  Because 
these  pressures  must  vary  like  e^s+y3K0  cosip)  on  ancj  Ko  k  u)/U0  in  very  low  Mach 
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number  flows,  the  incident  pressure  must  be  proportional  to  eMyn  very  close  to  S,  where  yn 
is  measured  outward  in  the  normal  direction  from  S  (i.e.  pj  must  decay  with  distance  ‘into’ 
S).  Therefore  (because  p0dvn/dt  =  -dpi/dyn  on  S),  we  can  write 


vn(k,K0  cos  ?/>,  co)  = 


-i\k\pi(k,  K0cosip,u)  _  -i\k\ps(k,  k0  cosip, u) 


p0u>  2  p0oj 

This  relation  permits  the  upwash  velocity  to  be  expressed  locally  in  terms  of  measured 
values  (or  an  empirical  model)  of  the  blocked  surface  pressure. 


(13) 


When  this  is  inserted  into  equation  (12),  and  the  result  is  summed  over  all  frequencies  u>, 
the  far  field  acoustic  pressure  is  obtained  in  the  form 


p(x,  t) 


sin  2  ^  sin(^/2)  ff°°  v/o;sgn(fc)^(fc, /t0cos^,o;)e  *^1+4 ) 


2v/2c^|x| 


//_ 


1*1* 


dkdu >, 


oo. 


(14) 


3.2  Curvature  noise  frequency  spectrum 

Let  a  finite  section  —\L<x3<\L  of  the  airfoil  span,  of  width  L  6  =  thickness  of 
the  turbulent  boundary  layer  at  O,  be  wetted  by  the  turbulent  flow.  When  the  flow  is 
statistically  stationary  in  time  the  blocked  surface  pressure  ps  satisfies 

(ps(k,  k3,(j)p*(k',  k3,u'))  «  -  u')5(k  -  k')P(k,k3,u),  L^>5,  (15) 

Z7T 

where  the  angle  brackets  (  )  represent  an  ensemble  average,  the  asterisk  denotes  complex 
conjugate,  and  P(k,  o>)  is  the  blocked  surface  pressure  wavenumber-frequency  spectrum  [17, 
21]. 

Then  if  the  one-sided  acoustic  pressure  frequency  spectrum  of  the  curvature  noise  is 
denoted  by  $o(x, oj),  such  that 

r  oo 

(p2(x,t))=  $0  (x,w)dw,  |x|->00,  (16) 

Jo 

equations  (14)  and  (15)  supply 


$o(x,w) 


wL  sin2  (6/ 2)  sin  ip  r°°  e 


87TC0|x 


i: 


i*i 


-P (k,  k0cos  ip,uj)  dk, 


oo. 


(17)- 


To  simplify  this  result  we  suppose  that  P (k,  k3,cu)  is  similar  in  form  to  the  wall-pressure 
spectrum  for  a  flat  wall  boundary  layer  [17,  21].  Thus,  if  U0  is  the  mean  flow  speed  parallel 
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to  S  just  outside  the  boundary  layer  at  0,  P (k,  fa,u ;)  will  have  a  very  large  maximum  in  the 
vicinity  of  k  =  u>/Uoc,  fa  =  0,  where  Uoc  ~  0.7 U0  is  a  boundary  layer  turbulence  convection 
velocity.  Now  k05  <  1,  because  the  Mach  number  of  the  flow  is  infinitesimal;  the  integral 
in  (17)  can  therefore  be  approximated  by  setting  k0  =  0  in  the  argument  of  P (k,  kq  cost/>,w) 
and  replacing  the  factor  e~lfc^/|fc|  by  its  value  at  k  =  u/Uoc,  i.e. 

roo  0 — | /c  | SR  £  q—oJ^/Uoc 

/  T-P(fc,  K0  cos  ip,  lo)  dk  &  — $PPM  —— ,  u  >  0,  (18) 

7—00  \k\  7T  \LOfUoc) 


where  4  U0/lo  is  the  Corcos  [22]  spanwise  correlation  length,  and  $pp(u;)  is  the  frequency 
spectrum  of  the  blocked  wall  pressure.  If  we  adopt  Chase’s  [21]  collation  of  experimental 
data,  we  can  model  the  wall  pressure  frequency  spectrum  by 


__  M.M)2 


{Poviy 


[(u><5*/C/0)2  +  <*p]5 


olv  =  0.12, 


(19) 


where  v%  is  the  boundary  layer  friction  velocity  at  O  [23]. 

Thus,  the  curvature  noise  spectrum  (17)  may  finally  be  cast  in  the  form 

(U/a)$ o(x,w)  _  0.7  /<U2  fon) 

(p0U2)2(La/\x\‘2)Msm2(d/2)smiP~  8ir2\a)  \U J  [(lj5*/U0)2  +  a2]§  ’  1 

where  M  =  U/c0,  and  where  v*  ~  0.035f/o,  <5*  are  the  local  values  of  the  friction  velocity 
and  boundary  layer  displacement  thickness  at  O. 
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4.  SLOT  NOISE 

4.1  Integral  representation  of  the  sound 

Sound  is  generated  at  the  slot  by  the  ‘scattering’  by  the  slot-lip  of  the  evanescent  (or 
hydrodynamic)  pressure  fluctuations  produced  by  turbulence  in  the  exterior  mean  flow 
boundary  layer  on  S  and  by  any  turbulence  within  the  high  speed  Coanda  jet.  To  estimate 
the  corresponding  contributions  to  the  far  field  acoustic  spectrum  the  thickness  of  the  slot 
lip  will  be  assumed  to  be  constant  and  equal  to  A,  and  the  influence  of  wall  curvature  in 
the  immediate  neighborhood  of  the  slot  will  be  ignored;  the  latter  assumption  is  justified  by 
the  observation  that  the  length  scales  of  the  dominant  slot  sources  must  be  comparable  to 
the  slot  width  h  which,  being  typically  smaller  than  about  10-2a,  is  itself  very  much  smaller 
than  the  local  radius  of  curvature  of  S.  This  simplified  analytical  model  of  the  Coanda  slot 
geometry  is  shown  in  Figure  4a;  for  the  purpose  of  calculating  the  interaction  of  the  flow 
with  the  lip  CD,  the  slot-duct  width  is  taken  to  be  uniform  and  equal  to  h.  In  the  figure 
the  ‘upper’  side  ED  represents  the  surface  S  of  the  hydrofoil  just  upstream  of  the  slot;  the 
end  A  of  the  lower  wall  B_A  is  the  continuation  of  S  downstream  of  the  slot;  the  upstream 
interior  ends  B±  of  the  duct  would  eventually  diverge  into  the  jet  plenum,  but  the  precise 
shape  of  the  duct  several  duct  widths  h  upstream  of  C  does  not  materially  affect  predictions 
of  the  sound. 

The  sound  produced  by  the  interaction  of  the  exterior  boundary  layer  and  the  Coanda 
jet  turbulence  with  the  lip  can  be  evaluated  from  the  general  formula  (7),  with  the  surface 
integral  confined  to  the  exterior  and  interior  surfaces  EDCB+  of  the  slot  lip,  and  to  the  wall 
B_A.  To  do  this  we  must  determine  the  form  of  the  potential  function  $*(y)  near  the  slot. 
This  function  satisfies  $*(y)  ~  $*(y)  (given  by  equation  (11))  at  distances  from  the  slot 
greatly  exceeding  the  slot  width  h.  Suppose  the  outer  corner  D  of  the  lip  (which  is  parallel 
to  the  span  of  the  hydrofoil)  is  at  Xs  =  {Xsi,Xs2),  then  near  the  slot  we  can  write 

(y)  =  $:(XS)  -  $oV*(y))  where  <  =  |V*;(X,)|.  (21) 

In  this  formula  represents  the  magnitude  of  the  potential  flow  velocity  at  Xs  defined 
by  the  velocity  potential  $*(y)  when  the  presence  of  the  slot  is  ignored.  Very  close  to  the 
hydrofoil,  this  velocity  is  parallel  to  the  nominal  surface  of  the  hydrofoil,  so  that  in  (21) 
the  function  ip*( y)  must  represent  the  distortion  of  this  hypothetical  flow  by  the  presence 
of  the  slot.  The  behavior  of  <p*( y)  can  be  specified  precisely  in  terms  of  local  coordinates 
(x,  y)  in  the  streamwise  xio^-plane,  with  origin  at  O  in  Figure  4a,  the  positive  rr-axis  lying 
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along  OA,  and  with  y  measured  outwards  into  the  fluid  from  0,  such  that  the  inner  and 
outer  corners  of  the  lip  C  and  D  are  respectively  at  (0,  h)  and  (0,  h  +  A).  In  other  words, 
in  the  vicinity  of  the  slot  we  use  the  local  coordinate  system  y  =  (x,y,y3).  Then  in  the 
exterior  region,  upstream  and  downstream  of  the  slot  in  Figure' 4a,  the  velocity  V<^*(y) 
must  represent  flow  at  unit  speed  in  the  positive  ^-direction. 

The  functional  form  of  <p*{ y)  =  <p*(x ,  y)  is  found  by  setting  Z  =  x  +  iy,  and  mapping  the 
upper  half  of  the  complex  Z-plane  exterior  to  the  semi-infinite  rectangular  domain  EDCB+ 
(i.e.  exterior  to  the  slot  lip  of  Figure  4a)  onto  the  upper  half  of  the  C-plane  by  means  of  the 
Schwarz-Christoffel  transformation  [24] 

t-  =  — ;=((C  +  1)5(C  +  h)*  +  (1  +  ln  [((  +  l)5  +  (C  + 11)2 

-  2^1n  [y/J(C  +  1)2  +(c  +  /^]  +\/juln<-  (1  -  v^)2ln  \A  -  > 

0  <  p  <  1,  Im  C  >  0.  (22) 

The  points  E,  D,  C,  B±,  A  all  map  onto  the  real  axis  in  the  C-plane  in  the  manner  indicated 
in  Figure  4b:  the  lip  EDCB+  maps  into  the  negative  real  axis  (— oo,  0),  and  lower  wall  B_A 
maps  onto  the  positive  real  axis.  The  value  of  p  is  given  in  terms  of  the  thickness  A  of  the 
lip  by 


The  image  in  the  C-plane  of  the  flow  past  the  slot  represented  by  the  potential  function 
<p*(y)  =  <p*{x ,  y)  is  a  flow  at  uniform  speed  in  the  positive  direction  of  the  real  C-axis,  whose 
complex  potential  is  therefore  a  linear  function  of  (.  According  to  (22),  Z  ~  at 

large  distances  from  the  slot  in  the  exterior  region,  so  that  unit  speed  at  large  distances 
from  the  slot  in  the  Z-plane  is  obtained  by  taking 

<p’{x>y)  =  Re  [— (C  +  1)1 .  (24) 

where  the  constant  T’  is  included  to  ensure  that  ip*  vanishes  at  the  corner  D  of  the 
lip.  By  evaluating  =  |V$*(X5)|  from  the  definition  (11),  we  may  now  cast  the  local 
representation  (21)  in  the  form 

r(y)  =  $-(Xl)  +  aRe(«Z)  +  l).  where  (25) 
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4.2  Passive  slot  noise  generated  by  the  external  boundary  layer 

‘Passive’  slot  noise  is  produced  by  the  interaction  of  the  exterior  mean  flow  boundary 
layer  turbulence  on  S  with  the  lip  of  the  jet  slot.  To  evaluate  the  surface  integral  of  equation 
(7)  for  this  component  of  the  sound  it  is  assumed  that  in  the  neighborhood  of  the  slot  lip 
all  of  the  turbulence  in  the  exterior  boundary  layer  lies  in  the  region  y  >  h  +  A  above  the 
lip.  Then  in  the  region  y  <  h  +  A,  each  frequency  component  pi(ytu)  of  the  hydrodynamic 
pressure  pi(y,t)  that  would  be  generated  by  the  boundary  layer  turbulence  when  diffraction 
at  the  surface  S  is  ignored  can  be  cast  in  the  form 

Pi(y,w)  =  /£>,  y  <  h  +  A,  (26) 

where  the  exponential  factor  e-(h+A~y)Vk2+kt  accounts  for  the  ‘evanescent’  decrease  in 
amplitude  of  the  induced  pressure  with  distance  h  +  A  —  y  ‘below’  the  boundary  layer. 
According  to  this  definition,  ps(k,k3,u )  =  2pi(k,  k3,u)  may  be  interpreted  as  the  blocked 
pressure  on  S  upstream  of  the  outer  corner  D  of  the  lip.  The  Fourier  transform  v(y,u>)  of 
the  upwash  velocity  on  S  may  now  be  calculated  from 

v(y,w)  =  Vpi(y,u>).  (27) 

ip0V 

Now  dS( y)  =  dsdy3  in  the  integral  fs  $*(y)un(y,  w)ei'Co|x_i3!/3|dS,(y)  of  equation  (7),  where 
s  is  the  curvilinear  distance  on  S  in  the  rcy-plane.  By  writing  Z  =  x  +  iy  and  noting  that 

f  ei{Ko\x-im\+km} dy3  =  2n6(k3  —  k0  cos  ip)em°^ ,  |x|  — >  oo, 

J  —  OO 

we  find 

f  .i.i  /y7rp^o|x|  TOO 

{  $*(y)vn{y,co)eiK°\x~'mldS(y)  = -  /  sgn(k) p s (k,  k0  cos  tp,u))ls(k)dk,  (28) 

J  S  Po^  •'-oo 

where 

( _ike-(h+A)\k\  r  ) 

ls(k)  =  C.C.  - f^*(y)e-%kzdz\,  k>  0, 

=  c.c.  {Zs(-A:)}  ,  k  <  0,  (29) 

and  the  notation  c.c.  {•}  implies  the  complex  conjugate  of  the  quantity  in  the  brace  brackets. 
The  integral  defining  Is  is  taken  over  the  piecewise  rectilinear  contour  7  =  AB_B+CDE 
in  the  Z-plane  of  Figure  4a,  starting  at  A  and  ending  at  E;  and  we  have  made  the  low 
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Mach  number  approximation  e  (h+A)^/k2+ni cos2  v  ~  e  IfcK/i+A)^  is  applicable  in  the 

hydrodynamic  domain  k  ~  u>/Us  »  u/c0  occupied  by  the  boundary  layer  pressure  sources. 

Only  those  contributions  to  the  integral  X s  from  the  neighborhood  of  the  mouth  of  the 
slot  are  relevant,  since  turbulence  sources  may  be  assumed  to  be  absent  very  far  upstream, 
and  to  die  out  (i.e.  become  incoherent  with  the  sources  near  the  slot)  at  large  distances 
downstream.  This  means,  that  when  using  the  approximation  (25)  for  $*(y)  (which  is 
valid  only  close  to  the  slot)  to  evaluate  Is,  it  is  permissible  to  secure  convergence  of  the 
integrals  by  integrating  by  parts  and  formally  discarding  contributions  from  Z  =  +00  and 
Z  =  -00  +  i(h  +  A).  By  noting  that  (  is  real  on  7,  we  then  find,  integrating  by  parts  once, 

Xs(k)  =  c.c.  |-e-(,l+A)|fc|  £  ^e"ifczdz|  =  c.c.  {e“(/l+A)|fc|  J°°  e-ifcZ(c)dc}  ,  k  >  0,  (30) 

The  final  integration  with  respect  to  (  is  taken  along  the  real  axis,  passing  just  above  the 
integrable  singularities  of  Z(Q. 

Hence,  collecting  together  these  results  for  (28)  and  substituting  into  equation  (7)  we 
finally  obtain,  after  summing  over  all  frequency  components,  the  passive  slot  noise  acoustic 
pressure: 


p(x,  t) 


— hsin2  ■0sin(0/2) 


(2tt)!  ^jIc-0\x\[(Xsl  +  3?/2)2  +  Xjy* 
x  Is(k)e~w^dkdu>,  |x|  — >  00. 


r  [[  Viujsgn(k)ps(k,K0cosip,u)) 

A  J  J  —  OO 


(31) 


This  can  be  used  to  calculate  the  corresponding  one-sided  acoustic  pressure  spectrum 
<hs(x>kO  of  the  passive  slot  noise  (defined  as  in  (16)): 


<Mx,oO 


coLh2  sin2  (9/2)  sin  ip 


8ttVo[(Xs1+3?/2)2  +  X22]5|x|2 


/OO 

| X${k)  |2P (&,  k0  cos  u)  dk 

-OO 


X 


OO, 


(32) 


where  P(k,  kz,u))  is  the  exterior  boundary  layer  blocked  surface  pressure  wavenumber- 
frequency  spectrum  just  upstream  of  the  slot.  When  kh  1  this  formula  reduces  in 
form  to  that  describing  the  production  of  sound  by  turbulent  boundary  layer  flow  over  a 
‘rectangular’  trailing  edge  of  thickness  A  when  the  flow  is  assumed  to  separate  at  the  corner 
adjacent  to  the  boundary  layer  [17].  In  this  limit  the  interaction  of  the  exterior  boundary 
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layer  turbulence  with  the  ‘image’  of  the  lip  in  the  inner  wall  B_A  of  figure  4a  becomes 
negligible. 

Numerical  estimates  can  be  made  by  evaluating  the  remaining  integral  in  (32)  by  an 
approximation  of  the  type  used  in  equation  (18),  together  with  a  corresponding  formula  (19) 
for  the  local  blocked  surface  pressure  frequency  spectrum  on  the  hydrofoil  just  upstream  of 
the  slot.  The  result  can  be  expressed  in  the  following  form  analogous  to  the  curvature  noise 
spectrum  (20) 

_ (E//q)$s(x,^) _ 

(po[/2)2(La/|x|2)Msin2(0/2)sint/> 

„  _1 _ ft  (hS.\  /M 4  2  (33) 

~  8tt5  [(xsl  +  5R/2)2  +  x22]l  V  a2  ;  VC/j  [M,/t/s)2  +  a2]§  ’ 

where  Usc  ~  0.7 Us,  Us  being  the  mean  flow  speed  indicated  in  Figure  4a  just  outside  the 
exterior  boundary  layer  above  the  slot,  and  6 *  and  u*  are  local  values  of  the  boundary  layer 
displacement  thickness  and  friction  velocity. 

In  general  \1s(<jv/USc)\2  must  be  evaluated  numerically  from  the  definition  (30).  Because 
the  integral  in  (30)  is  dominated  by  contributions  from  the  neighborhood  of  the  lip, 
numerical  evaluation  is  facilitated  (for  k  >  0)  by  first  deforming  the  contour  in  the  £-plane 
onto  the  piecewise  rectilinear  path  joining  the  points  £  =  — 1  —  *oo,  —1,  0,  -zoo,  where  the 
section  along  the  real  axis  between  £  =  —1,  0  passes  above  the  singularities  of  Z( £).  On 
this  contour  the  integrand  is  exponentially  small  at  the  two  extremes  where  £  -4  —zoo.  The 
following  two  limiting  cases  can  be  evaluated  explicitly: 


\Uk)\2=l 


|A;|A  »  1,  A^0; 
A  =  0. 


(34) 


These  formulae  respectively  define  the  high  frequency  behavior  when  the  lip  thickness 
A  ^  0,  and  the  exact  behavior  for  all  frequencies  when  A  =  0  (i.e.  when  A  </i). 

4.3  Slot  noise  generated  by  the  Coanda  jet 

In  typical  underwater  applications,  the  Coanda  jet  Reynolds  number  Ujh/v  ~  105 
{y  ~  10-6  m2/sec  being  the  kinematic  viscosity  of  water).  This  would  normally  imply 
that  the  jet  flow  from  the  slot  is  fully  turbulent.  However,  the  exit  flow  from  the  plenum 
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converges  towards  the  slot  thereby  inhibiting  instability  and  the  onset  of  turbulence,  so  that 
some  care  must  be  exercised  in  estimating  the  turbulence  properties  of  the  jet.  In  undersea 
applications  fouling  of  interior  surfaces  may  be  a  source  of  turbulence. 

By  symmetry  it  may  be  assumed  that  the  blocked  surface  pressures  ps  produced  by 
the  turbulent  jet  on  the  upper  and  lower  slot  duct  walls  are  statistically  the  same.  Then, 
to  evaluate  the  surface  integral  §s  $*(y)vn(y,  (v)etKo^x~i3y^dS(y)  of  equation  (7)  by  the 
procedure  described  above,  we  can  take  the  ‘incident’  hydrodynamic  pressure  Pi(y,w)  to  be 
given  by 


Pi(y,w) 


//_> 

//_> 


,  Jfe3,  u)Skx+k*y*)-{y-h)'/^dkdk 3,  y  >  h  on  B+CDE, 
,  A&,  u)ei{kx+k3y3)+y^+^dkdk3,  y<  0  on  B_A, 


(35) 


where  in  these  two  formulae,  pi(k,k3,u>)  respectively  represents  the  Fourier  transform  of 
the  incident  pressure  generated  by  the  jet  turbulence  on  B+C  and  B_A.  In  each  case  the 
upwash  velocity  is  calculated  from  (27).  The  radiated  pressures  may  now  be  calculated 
by  the  method  used  in  §4.2  for  the  passive  slot  noise.  In  doing  this,  however,  it  should 
be  noted  that  the  wall  pressures  on  the  upper  and  lower  slot  duct  walls  will  normally  be 
statistically  independent  except  at  very  low  frequencies.  The  net  acoustic  pressure  spectrum 
will  therefore  be  equal  to  the  sum  of  the  separate  spectra  generated  by  the  upper  and  lower 
Coanda  jet  boundary  layers. 


In  the  acoustic  far  field  we  find 

±h  sin ^  ip  sin(6/2) 


P±(x,t) 


(2tt)2  ^jjT0\x\[(Xsl  +  3ty2)2  +  X2s2] 


1 

4 


x  f  Vivsgn(k)ps(k,K0cosip,ui)l±(k)e  lu^dkdu,  |x|  — >  oo,  (36) 

where  p±(x.,  t)  denotes  the  acoustic  pressure  respectively  attributable  to  the  upper/lower 
boundary  layers  of  the  Coanda  jet,  and 

X+{h)  =  <fE  $*(y )eikzdZ  =  elfc|/l  [°  eifcz(c)dC,  *  >  0, 

OL  Jb+  J— oo 

=  c.c.  {I+(—k)} ,  A:  <  0,  (37) 

I L(Jfc)  =  c.c.|^£"r(y)e-ifc^|  =  c.c.  { ^°°  e"ifcZ(c)dc}  ,  ^>0, 

=  c.c.  {!_(—&)} ,  k  <0,  (38) 


17 


Report  No.  AM  00-004 


Boston  University,  College  of  Engineering 


The  corresponding  one-sided  acoustic  pressure  spectra  <3>±(x,a;)  are  given  by 

\T±(k)\2P(k,K0  cos  ip,  oj)  dk,  |x|  — >  oo, 

(39) 

where  P(k,k3,cj )  is  the  blocked  pressure  wavenumber-frequency  spectrum  for  the  jet, 
assumed  to  be  the  same  on  the  upper  and  lower  walls  of  the  jet  slot.  The  remaining 
wavenumber  integral  can  be  approximated  as  previously,  by  expanding  about  the  peak  of 
P(k, K0costp,uj)  «  P(/c,0,u;)  at  k  =  u/UjC  (where  UjC  ~  0.7Uj).  Using  the  corresponding 
representation  (19)  of  the  blocked  pressure  frequency  spectrum,  we  then  find 


ujLh?  sin2  (9/ 2)  sin  ip  r°° 

8tt^c0[(Xs1  +  R/2)2  +  Xiy5|x|2  J-<*> 


_ (U/q)^(x^) _ 

(p0U2)2(La/\x.\2)M  sin2(9/2)  sin  ip 

„  J _ h  ( hS A  (v*\ 4  W (40' 

~  8^5  [(*.i  +  ft/2)2  +  X22)\  \a2)\u)  [(uS./Uj)2  +  og]f  ’ 

where  8 *  and  u*  are  the  boundary  layer  displacement  thickness  and  friction  velocity  in  the 
Coanda  jet  flow  from  the  slot. 

The  values  of  \l±(co /Ujc)\2  must  be  found  by  numerical  integration,  as  before.  However, 
the  following  special  cases  may  be  noted: 


\X+{k)\2  = 


3F2  (|) /  J 3  >  WA»1>  A^°; 

2'n2  j  \k\h  (l  —  e-2lfcl/l)] ,  A  =  0; 
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5.  NUMERICAL  RESULTS 

To  illustrate  predictions  of  these  formulae  we  consider  a  typical  application  involving  a 
nominally  elliptic  hydrofoil  at  zero  mean  angle  of  attack  with  semi-major  and  minor  axes 
a  =  1.5  m  and  b  =  0.2a  respectively,  and  with  a  rounded  trailing  edge  of  radius  9?.  The 
following  parameter  values  are  representative  of  what  might  be  encountered  when  a  Coanda 
jet  is  used  to  generate  a  lift  coefficient  Cl  ~  3.5  in  water  at  a  forward  speed  of  U  =  5  knots 
(2.5  m/sec): 

h  =  0.006a 
ft  =  0.08a 

{(Xsl  +  ft/2)2  +  X22}5  =  0.1a 

U0  =  2.8  U 

Uj  =  6  U 

Us  =  OAUj 

It  is  assumed  that  entrainment  by  the  jet  produces  a  local  increase  in  the  exterior  mean 
velocity  Us  near  the  slot;  the  choice  Us  =  OAUj  is  based  on  the  experimental  data  of  Novak 
et  al.  [8  -  10].  For  the  exterior  boundary  layer  flow  we  can  estimate  that  near  the  slot  and 
the  rounded  trailing  edge  the  displacement  thickness  5 *  «  0.004a.  Under  ideal  conditions 
the  Coanda  jet  is  laminar  just  upstream  of  the  slot  exit;  we  call  this  case  (a),  in  which 
there  is  no  direct  contribution  to  the  sound  from  the  interaction  of  the  jet  with  the  slot 
lip.  However,  according  to  the  observations  of  Novak  et  al.  [8  -  10]  turbulence  velocities 
of  magnitude  ~  0.1  Uj  occur  near  the  slot  exit,  where  the  same  data  indicates  that  the 
boundary  layer  thickness  ~  0.15 h\  we  call  this  case  (b)  and  take  6 *  «  0.15/i/8  =  0.019/i 
[23].  The  friction  velocity  is  taken  to  be  related  to  the  local  mean  stream  velocity  V  just 
outside  the  corresponding  turbulent  boundary  layer  by  v *  =  0.035V  (we  take  V  =  Uj  in  the 
Coanda  jet)  [23]. 

These  values  and  equations  (20),  (33)  and  (40)  may  now  be  used  to  calculate  the  acoustic 
pressure  frequency  spectrum  for  cases  (a)  and  (b)  for  prescribed  values  of  the  lip  thickness 
A.  Thus,  the  solid  curves  in  Figures  5  and  6  represent  predictions  of  the  overall  sound 
pressure  spectrum 

(poU2)2(La/|x|2)Msin2(0/2)  sin^» 

in  these  two  cases  for  a  lip  of  zero  thickness  (A  =  0)  where,  respectively  for  cases  (a)  and  (b), 
4>(x,u;)  =  $o(x,w)  +  $s(x,  w)  and  <&(x,cj)  =  4>0(x) oj)  +  $s(x)w)  +  (x,w)  +  ^(xjw). 
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The  broken  line  curves  in  the  figure  represent  the  individual  contributions  from  the 
curvature  noise  $0,  the  passive  slot  noise  $s  and  (in  Figure  6)  the  Coanda  jet-slot 
interaction  <E>j\  The  very  low  frequency  sound  is  dominated  by  the  curvature  noise.  At 
high  frequencies  the  passive  slot  noise  is  dominant  in  case  (a),  but  Figure  6  shows  that 
when  turbulence  is  present  in  the  Coanda  jet  the  spectrum  is  extended  out  to  much  higher 
frequencies  attributable  to  the  interaction  of  the  jet  turbulence  with  the  lip.  In  case  (b), 
and  for  these  particular  parameter  values,  the  relative  contribution  from  the  passive  slot 
noise  is  negligible.  In  both  cases  $(x,u>)  is  predicted  to  decrease  with  /  =  uj/2-n  like 
<&pp(uj)/lj  ~  /~2  at  high  frequencies. 

At  frequencies  below  about  500  Hz  the  assumption  used  in  the  calculation  that  the 
acoustic  wavelength  is  large  compared  to  the  hydrofoil  chord  ceases  to  be  valid,  and  results 
at  these  frequencies  should  strictly  be  corrected  by  application  of  the  interpolation  factor 
given  in  [20].  The  correction  would  decrease  the  amplitude  of  the  sound  predicted  at  lower 
frequencies,  because  of  the  reduced  efficiency  with  which  sound  is  produced  by  a  hydrofoil 
of  compact  chord.  The  magnitude  of  the  low  frequency  correction  depends  strongly  on  the 
radiation  direction  but,  because  the  fluid-surface  interaction  mechanism  is  independent  of 
the  chord  for  small  scale  turbulence,  we  can  obtain  an  overall  picture  of  the  importance 
of  trailing  edge  rounding  and  the  slot-generated  sound  at  these  frequencies  by  comparing 
our  present  predictions  for  a  non-compact  chord  with  the  corresponding  predictions  for  a 
flat-plate  hydrofoil  of  the  same  chord,  with  the  same  turbulent  flow  over  its  trailing  edge. 
For  such  a  flat-plate  hydrofoil  the  mean  flow  velocity  outside  the  boundary  layer  is  constant 
and  equal  to  U,  and  for  turbulent  flow  on  one  side  only,  we  find  [17] 

_ (U/a)${x.,u) _ „  Q-7  (S*Y  (VA*  {ut*/U)  (u) 

(po[/2)2(La/|x|2)Msin2(0/2)sin^  ~  2tt2  \a  )  \U )  [(cjSJU)2  +  a2]§ '  V  ' 

This  is  plotted  as  the  dotted  curves  in  Figures  5  and  6.  Comparison  with  the  spectrum 
for  the  overall  CC-hydrofoil  noise  reveals  a  broad  level  of  agreement  at  very  low  frequencies, 
and  also  at  high  frequencies  when  the  Coanda  jet  is  turbulent-free.  However,  at  high 
frequencies  the  radiation  generated  by  the  turbulent  Coanda  jet  is  typically  20  dB  in  excess 
of  the  flat-plate  edge  noise. 

Figures  7  and  8  display  the  corresponding  predictions  respectively  for  cases  (a)  and  (b) 
when  the  slot  lip  has  thickness  A  =  0.5 h.  In  case  (a)  (no  jet  turbulence)  the  passive  slot 
noise  is  smaller  (relative  to  the  case  A  =  0)  by  about  5  dB  at  low  frequencies;  the  reduction 
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is  larger  at  high  frequencies  because  the  finite  thickness  of  the  lip  causes  the  spectrum 
to  decrease  like  §pp{w)/u)3  ~  /“I  (instead  of  like  /“ 2  when  A  =  0);  the  efficiency  of 
production  of  sound  at  the  slot  is  in  this  case  smaller  than  for  the  corresponding  flat-plate 
edge  flow  represented  by  the  dotted  curve.  When  the  Coanda  jet  is  turbulent  (Figure  8, 
Case(b))  the  higher  frequencies  are  once  again  dominated  by  the  jet-slot  interaction  noise 
(which  also  decreases  like  /~s). 
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6.  CONCLUSION 

Coanda  wall  jet  circulation  control  permits  the  lift  generated  by  a  hydrofoil  to  be  greatly 
increased  at  low  mean  flow  speeds  without  the  deployment  of  mechanical  lift  augmentation 
devices  such  as  a  trailing  edge  flap.  However,  because  of  the  relatively  high  jet  speeds 
required  to  maintain  attached  flow  around  the  hydrofoil,  additional  self-noise  is  produced 
over  and  above  that  usually  attributed  to  the  boundary  layer  turbulence  (generated  by 
boundary  layer  instability  even  when  the  incident  mean  flow  is  nominally  steady)  interacting 
with  the  trailing  edge. 

Three  noise  sources  that  are  likely  to  dominate  the  production  of  sound  in  the  vicinity 
of  the  trailing  edge  of  the  hydrofoil  have  been  identified.  Theoretical  predictions  of  their 
relative  contributions  to  the  overall  self-noise  of  the  CC-hydrofoil  are  critically  dependent 
on  a  proper  understanding  of  the  turbulence  characteristics  of  the  flow,  especially  within 
the  Coanda  jet.  The  following  quantitative  conclusions  are  based  on  particular  numerical 
estimates  for  these  quantities  used  in  this  report,  and  should  therefore  be  regarded  as 
tentative. 

Very  low  frequencies  are  dominated  by  ‘curvature  noise’,  generated  by  the  interaction  of 
boundary  layer  turbulence  with  the  rounded  trailing  edge  of  the  CC-hydrofoil;  it  is  similar  in 
character  and  amplitude  to  the  low  frequency  component  of  the  trailing  edge  noise  produced 
by  a  hydrofoil  of  equal  chord,  but  with  a  sharp  trailing  edge.  The  jet  slot  is  the  main 
source  at  higher  frequencies.  ‘Passive  slot  noise’  is  produced  by  the  ‘scattering’  into  sound 
of  the  near  field  pressure  fluctuations  of  turbulence  in  the  exterior  mean  boundary  layer 
flow  over  the  slot.  This  is  much  stronger  than  the  curvature  noise  at  higher  frequencies, 
and  is  comparable  in  magnitude  to  the  sharp  edged  trailing  edge  noise.  However,  if  the 
Coanda  jet  is  turbulent,  the  sound  produced  by  its  interaction  with  the  lip  of  the  slot  tends 
to  be  the  dominant  source  at  high  frequencies  (say  for  fh/Uj  >  1).  The  overall  level  of 
the  interaction  noise  depends  on  the  thickness  A  of  the  slot  lip:  when  A  =  0  the  Coanda 
jet  noise  at  higher  frequencies  is  typically  20  dB  or  more  greater  than  conventional  trailing 
edge  noise.  The  difference  becomes  smaller  as  the  lip  increases  in  thickness,  since  this  tends 
to  reduce  the  efficiency  of  the  edge  scattering. 

It  has  been  implicitly  assumed  that  the  acoustic  impedance  of  the  jet  slot  is  infinitely 
large.  This  means  that  no  account  is  taken  of  ‘cavity  resonances’  of  the  jet  plenum. 

Such  resonances  probably  occur  in  practice  at  very  low  frequencies.  In  principle  their 
contributions  can  be  included  by  modifying  the  hydroacoustic  Green’s  function  in  the 
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manner  described  in  [25],  although  the  modification  would  depend  on  a  detailed  specification 
of  the  plenum  characteristics  and  would  be  useful  only  if  observation  indicated  the  presence 
of  resonance  peaks  in  the  acoustic  spectrum. 


23 


Report  No.  AM  00-004 


Boston  University,  College  of  Engineering 


REFERENCES 

1.  H.  Coanda  1932  Brevet  Invent.  Gr.  Cl.  2  No.  762688  Republique  Frangaise.  Procede 
de  propulsion  dans  un  fluide. 

2.  P.  W.  Carpenter  and  P.  N.  Green  1997  Journal  of  Sound  and  Vibration  208,  777 
-  802.  The  aeroacoustics  and  aerodynamics  of  high-speed  Coanda  devices,  Part  I: 
Conventional  arrangement  of  exit  nozzle  and  surface. 

3.  L.  J.  Pritchard  1957  Journal  of  the  Royal  Aeronautical  Society  61,  149  -  180.  The 
dawn  of  aerodynamics. 

4.  R.  Wille  and  H.  Fernholz  1965  Journal  of  Fluid  Mechanics  23,  801  -  819.  Report  on 
the  first  European  Mechanics  Colloquium  on  the  Coanda  effect. 

5.  R.  J.  Englar  and  C.  A.  Applegate  1984  Circulation  control  -  A  bibliography  of  DTNSRDC 
research  and  selected  outside  references  (Jan  1969  -  Dec.  1983).  David  Taylor  Naval 
Ship  Research  and  Development  Center  Report  No.  DTNSRDC  84/052,  Bethesda  MD. 

6.  I.  C.  Cheeseman  and  A.  R.  Seed  1966  Journal  of  the  Royal  Aeronautical  Society  71, 
451  -  467.  The  application  of  circulation  control  by  blowing  to  helicopter  rotors. 

7.  N.  J.  Wood  and  J.  N.  Nielsen  1986  American  Institute  of  Aeronautics  and  Astronautics 
Journal  23,  865  -  875.  Circulation  control  airfoils  as  applied  to  rotary-wing  aircraft. 

8.  Novak,  C.  J.  and  Cornelius,  K.  C.  1986  Lockheed-Georgia  Company  Report  L87R0618 
Volume  II:  Advanced  (circulation  control)  airfoil  performance  and  analysis:  experi¬ 
mental  program  summary. 

9.  Novak,  C.  J.  and  Cornelius,  K.  C.  1986  Lockheed-Georgia  Company  Report  L87R0618 
Volume  III:  Advanced  (circulation  control)  airfoil  performance  and  analysis:  experi¬ 
mental  program  data. 

10.  Novak,  C.  J.,  Cornelius,  K.  C.  and  Roads,  R.  K.  1987  AIAA  Paper  87-0155.  Experi¬ 
mental  investigation  of  the  circular  wall  jet  on  a  circulation  control  airfoil. 

11.  J.  N.  Nielsen  (editor)  1987  Proceedings  of  the  Circulation  Control  workshop  1986, 
NASA  Conference  Publication  2432. 

12.  A.  Zandieh  and  J.  G.  Leishman  1993  American  Institute  of  Aeronautics  and  Astronau¬ 
tics  Journal  31,  1769  -  1776.  Boundary  layer  and  pressure  measurements  on  a  circular 
cylinder  with  unsteady  circulation  control. 


24 


Report  No.  AM  00-004 


Boston  University,  College  of  Engineering 


13.  V.  Raghavan,  I.  Chopra  and  S.  Pai  1988  Journal  of  the  American  Helicopter  Society 
33,  28  -  37.  Circulation  control  airfoils  in  an  unsteady  flow. 

14.  M.  Salikudden,  W.  H.  Brown  and  K.  K.  Ahuja  1987  Journal  of  Aircraft  24,  55  -  64. 
Noise  from  a  circulation  control  wing  with  upper  surface  blowing. 

15.  P.  W.  Carpenter  and  C.  Smith  1997  Journal  of  Sound  and  Vibration  208,  803  -  822. 
The  aeroacoustics  and  aerodynamics  of  high-speed  Coanda  devices,  Part  II:  Effects  of 
modifications  for  flow  control  and  noise  reduction. 

16.  M.  J.  Lighthill  1952  Proceedings  of  the  Royal  Society  A211,  564  -  587.  On  sound 
generated  aerodynamically.  Part  I:  General  theory. 

17.  M.  S.  Howe  1998  Acoustics  of  Fluid- Structure  Interactions ,  Cambridge  University 
Press. 

18.  M.  S.  Howe  1999  Journal  of  Sound  and  Vibration  225,  211  -  238.  Trailing  edge  noise 
at  low  Mach  numbers. 

19.  G.  K.  Batchelor  1967  An  Introduction  to  Fluid  Dynamics.  Cambridge  University  Press. 

20.  M.  S.  Howe  1999  Quarterly  Journal  of  Mechanics  and  Applied  Mathematics  (in  press) 
(also,  Boston  University  Department  of  Aerospace  and  Mechanical  Engineering  Report 
No.  AM99-010).  Edge-source  acoustic  Green’s  function  for  an  airfoil  of  arbitrary  chord, 
with  application  to  trailing  edge  noise. 

21.  D.  M.  Chase  1980  Journal  of  Sound  and  Vibration  70,  29  -  67.  Modeling  the  wavevector- 
frequency  spectrum  of  turbulent  boundary  layer  wall  pressure. 

22.  G.  M.  Corcos  1964  Journal  of  Fluid  Mechanics  18,  353  -  378.  The  structure  of  the 
turbulent  pressure  field  in  boundary  layer  flows. 

23.  J.  O.  Hinze  1975  Turbulence  (second  edition).  New  York:  McGraw-Hill. 

24.  Horace  Lamb  1932  Hydrodynamics  (6th.  ed.).  Cambridge  University  Press  (Reprinted 
1993). 

25.  M.  S.  Howe  1976  Journal  of  Sound  and  Vibration  45,  427  -  440.  On  the  Helmholtz 
resonator. 


25 


Report  No.  AM  00-004 


Boston  University,  College  of  Engineering 


Figure  1.  Schematic  Coanda  jet  circulation  control  hydrofoil,  set  at  zero  angle  of  at¬ 
tack  to  a  nominally  uniform  mean  flow  at  speed  U ,  showing  the  jet  slot  and 
plenum. 
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Figure  2.  (a)  Trajectory  of  the  line  vortex  (5)  around  the  trailing  edge  of  an  elliptic  airfoil  of  semi- 
major  and  minor  axes  a,  b  with  b/a  =  0.2,  for  U  =  2.5  m/sec.,  a  =  1.5  m,  CL  =  3.5. 
The  vortex  is  released  at  the  midchord  position  at  a  standoff  distance  equal  to  the 
nominal  boundary  layer  displacement  thickness  (5*  =  0.0024a.  (b)  Nondimensional  far 
field  acoustic  pressure  p(x,  t)/[(p0U2)(FLcos  ©/c0a|x|)]  plotted  as  a  function  of  the 
retarded  time  [t]  =  t  —  |x|/c0,  showing  the  ‘curvature  noise’  produced  when  the  vortex 
passes  around  the  trailing  edge  at  O,  and  the  ‘separation  noise’  generated  near  A. 


27 


Report  No.  AM  00-004 


Figure  4.  (a)  Local  geometry  in  the  Z-plane  of  the  Coanda  jet  slot  and  lip,  with. 

Coanda  jet  velocity  Uj  and  exterior  mean  flow  speed  Us]  (b)  Image  in 
Im  C  >  0  of  the  fluid  region  (exterior  to  the  lip)  in  the  upper  half  of  the 
Z-  plane. 
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frequency  Hz 


Figure  5.  Overall  acoustic  pressure  frequency  spectrum  ( — 

(U/a)$(x,u) 


10  x  log10 


G 


(p0U2)2 (La /|x|2)M  sin2 (9/2)  simp) 

and  the  component  spectra  $o,  $s  ( - )  for  conditions  (43)  for  a 

turbulence-free  Coanda  jet  (Case  (a))  and  a  lip  of  thickness  A  =  0.  The 
dotted  curve  (•  •  •)  is  the  corresponding  edge  noise  spectrum  (44)  for  a 
flat- plate  hydrofoil  of  equal  chord. 
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Figure  6.  Overall  acoustic  pressure  frequency  spectrum  ( - ) 

_ (U/am^Qj) _ 

(po£/2)2(La/|x|2)Msin2(0/2)  simp 

and  the  component  spectra  $o,  $S)  ( - )  for  conditions  (43)  for  a 

turbulent  Coanda  jet  (Case  (b))  and  a  lip  of  thickness  A  =  0.  The  dotted 
curve  (•  •  •)  is  the  corresponding  edge  noise  spectrum  (44)  for  a  flat-plate 
hydrofoil  of  equal  chord. 
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Figure  7.  Overall  acoustic  pressure  frequency  spectrum  ( - ) 

_ (U/a)$(x,u) _ 

(, p0U2)2(La/\x\2)M  sin2  (0/2)  sin  ip 

and  the  component  spectra  4>o,  $s  ( - )  f°r  conditions  (43)  for  a 

turbulence-free  Coanda  jet  (Case  (a))  and  a  lip  of  thickness  A  =  0.5 h. 
The  dotted  curve  (•  •  •)  is  the  corresponding  edge  noise  spectrum  (44)  for 
a  flat-plate  hydrofoil  of  equal  chord. 
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Figure  8,  Overall  acoustic  pressure  frequency  spectrum  ( - ) 

(  (U/a)$(x,uj) 

10  X  log10  ^poC/2)2(Lo/|x|2)Msin2^/2)Sin^ 

and  the  component  spectra  $o,  $s,  <hj:  ( - )  f°r  conditions  (43)  for 

a  turbulent  Coanda  jet  (Case  (b))  and  a  lip  of  thickness  A  =  0.5 h.  The 
dotted  curve  (•  •  •)  is  the  corresponding  edge  noise  spectrum  (44)  for  a 
flat-plate  hydrofoil  of  equal  chord. 
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